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Abstract 

New problem is considered that is to find nonlinear difi^erential equa- 
tions with special solutions. Method is presented to construct nonlinear 
ordinary differential equations with exact solution. Crucial step to the 
method is the assumption that nonlinear differential equations have ex- 
act solution which is general solution of the simplest integrable equation. 
The Riccati equation is shown to be a building block to find a lot of non- 
linear differential equations with exact solutions. Nonlinear differential 
equations of the second, third and fourth order with special solutions are 
given. Most of these equations are used at the description of processes in 
physics and in theory of nonlinear waves. 
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PACS: 02.30.Hq - Ordinary differential equations 

1 Introduction 

In recent years one can observe a splash of papers where authors presented a 
lot of different approaches to look for exact solutions of nonlinear differential 
equations m21Eli4liSliiSn2llS»l2j • There are two reasons to make the study in 
this direction. First there is a great interest to the investigation of nonlinear 
phenomenon. Secondly we have codes Maple, Mathematica and other ones to 
conduct a lot of symbolical calculations. 

It is well known that all nonlinear differential equations can be connectionally 
divided into three types: exactly solvable, partially solvable and those that have 
no exact solution. 

Consider nonlinear evolution equation 

El[u] EE El{u,Uf,Ux, = (1.1) 

Assume we need to have exact solutions of this equation. First of all one can 
try to solve a question about the integrability of this equation. For this aim one 
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can apply the Painleve test to check the integrabihty of this equation JOI- As 
a result we can have two variants. First variant is the equation pass the 
Pailevee test and we have the necessary condition for integrabihty of equation 
in this case. Second variant is the equation does not pass the Painleve 

test. In this case we can look for some transformation to obtain equation that 
can be passed the Painleve test. However we often obtain that the equation 
does not pass the Painleve test and we have not got any transformation 
to obtain good form for the origin equation. Unfortunately a lot of nonlinear 
evolution equation does not pass the Painleve test. In any case one can search 
exact solutions of nonlinear differential equation. In fact one can look for exact 
solutions of nonlinear differential equation of all types without solution of the 
problem about the integrabihty. 

Usually we look for exact solution of nonlinear evolution equation taking 
into account the travelling wave and search exact solution of equation in 
the form 



u{x,t) = y{z), z ~ X — C'ot (1-2) 

As a result we have that the equation reduces to the nonlinear ordinary 
differential equation (ODE) 



E2[y] = E2{y,y.,...,z)^0 (1.3) 

To obtain exact solutions of equation 1)1. 3|1 one can apply different approaches 
[aOliniiniinilTllTnilTniiniliHlIiniEni. However one can note that the most 
methods that are used to search exact solutions take into account the singular 
analysis for solutions of the nonlinear differential equations |lUllllll2ll21l22l23j . 

Even though the investigators are not aware of the analytical theory of nonlinear 
differential equations they apply as a rule its ideas and approaches. 

Using the singular analysis first of all we have to consider the leading mem- 
bers of equation . After that we find the singularity for solution of equation 
(|1.3|) . Further the truncated expansion is used to have the transformation to 
search exact solutions of nonlinear ODEs. At this point one can use some trial 
functions (hyperbolic and elliptic and so on) to look for exact solutions of non- 
linear ODEs ^A.M2^.M.2LM,M.M,- 

However one can note that hyperbolic and elliptic functions are general so- 
lutions of nonlinear exactly solvable equations. We have as a rule that partially 
solvable nonlinear differential equations have exact solutions that are general 
solutions of solvable equations of lesser order. 

In this paper we are going to start the solution of the new problem that is 
to find nonlinear ordinary differential equations of polynomial form which have 
special solutions. Let us explain the idea of our work. It is well known there 
is the great problem that is to find nonlinear differential equations that are 
integrable. It is very important because we often want to have general solution 
of nonlinear ordinary equations. For a example Painleve and his school found 50 
canonical forms of the second order class for the ordinary differential equations 
that are integrable equations. However sometimes we can content ourselves with 
some special solutions because a lot of differential equations are nonintegrable 
ones although they are intensively used in physics and look as simple equations. 
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It is important to find some special solutions of these equations that are called 
exact solutions. In this connection it is convenient to have the list of nonlinear 
differential equations that are not integrable but have special solutions. 

The aim of this work is to present new method to find nonlinear differential 
equations with exact solutions. Using our approach we give a class of nonlinear 
ordinary differential equations that have exact solutions. Special solutions of 
our nonlinear ordinary differential equations are found via general solution of 
the Riccati equation. 

As this takes place we do not solve the problem of finding all possible exact 
solutions for our class of nonlinear ordinary differential equations. Our class 
contains a lot of exactly solvable equations. In this paper we give only the poly- 
nomial class of nonlinear ordinary differential equations with special solutions. 

The outline of this paper is as follows. In section 2 we present the method 
to find nonlinear ordinary differential equations (ODEs) with special solutions. 
Nonlinear ODEs with exact solutions of the first, second, third and fourth degree 
singularities are given in section 3, 4, 5 and 6. Example of nonlinear ODE with 
exact solution of the fifth degree singularity is considered in section 7. 

2 Method applied 

Let us discuss the method that can be applied to find nonlinear differential 
equations with exact solutions. One can note that most nonlinear ordinary 
differential equations has exact solutions that are general solutions of differential 
equations of lesser order. Much more than that exact equations for the most 
nonlinear differential equations are determined via general solution of the Riccati 
equation. This is so indeed because one can note that the most approaches to 
search exact solutions of nonlinear ordinary differential equations are based on 
general solution of the Riccati equation. The application of the tanh method 
confirms this idea PHlHlEniEIllIll ■ This is due to the fact that the Riccati 
equation is simplest equation of first order of polynomial form that have the 
Painleve property |TT)lll7ll21j . 

The Riccati equation can be presented in the form 

^ + + = (2.1) 
We have the following simple theorem. 
Theorem 2.1. Let Y{z) he solution of equation H2.1|l than equations 

= 2Y^ + 2qY - q, (2.2) 



Y,,, = -6Y^ - 8qY^ + 2q,Y - 2q^ - q,, (2.3) 

y,,,, - 24r-^ + 40qy^ + IQq^Y - lOq.Y^ + 2q,,Y - 6qq, - q,,, (2.4) 
have special solutions that are expressed via general solution of equation H2.1|l 
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Proof. Theorem 2.1 is proved by differentiation of (|2.1|) with respect to z and 
substitution Yz from equation (|2.1() and so on into expressions obtained. 

□ 

Corollary 2.1. Equations 

^ 2Y^ - 2aY (2.5) 
Y,,, = -6Y^ + 8aY^ - 2o? (2.6) 

Yy_zzz = 24F^ - 40ar^ + \%c?Y (2.7) 

have special solutions Y{z) where Y{z) is the general solutions of the Riccati 
equation in the form 

n + y2 _ = (2.8) 
It is well known that general solutions of equation (|2.8(l is 

y(z) = ^/atanh (Vaz + (/3o) (2.9) 

where is arbitrary constant. 

We want to find nonlinear ordinary differential equation which have special 
solutions that are determined via general solution of the Riccati equation. 

Algorithm of our method can be presented by four steps. At the first step 
we choose the singularity of special solution and give the form of this solution. 
At the second step we take the order of nonlinear ordinary differential equation 
what we want to search. The third step lies in the fact that we write the general 
form of nonlinear differential equation taking into account the singularity of 
the solution and the given order for nonlinear differential equation. The fourth 
step contains calculations. As a result we find limitations for the parameters 
in order for nonlinear differential equation has exact solutions. At this step we 
have nonlinear ODE with exact solutions. 

Let us demonstrate our approach. With this aim let us find a nonlinear 
ordinary differential equation of the second order with solution of the second 
degree singularity. This solution takes the form 

y{z) ^ Ao + A,Y + A^Y^ (2.10) 

Here Y{z) is a solution of equation (|2.8|) . Without loss of generality we can take 
A2 = -12. 

First of all let us write the general form of the nonlinear second order ordi- 
nary differential equation with solution H2.10|l . It takes the form 

El = y,, + aoy^ - boVz - Coy + Ci = (2.11) 
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Here oq, 5o, Cq and Ci are unknown parameters of equation 
This equation was written using singularity of solution (|2.10() . Actually the 
singularity of the first term in (|2.11l) is equal to 4 and for the polynomial form 
of the equation we have to take into account the same singularity of nonlin- 
ear expression for the second term (|2.11(l . Other terms in H2.11|l have lesser 
singularity. 

We want to find equation H2.11|l that has exact solution H2.9|l . We can do 
this if we define values of parameters ao,6o,Co and Ci. Substituting (|2.10() at 
A2 = —12 into equation (|2.1H) and taking into account equations l|2.5|l and H2.8|l 
we have equation in the form 

(-72 + 144 ao) + (-24 ao^i + 2 Ai - 24 60) 
+ (-24 aoAo + aoAi^ + 12 Co + boAi + 96 a) 

(2.12) 

+ (-C0A1 - 2 Aia + 2 ao^o^i + 24 60a) Y~ 

-24 + aoAo^ - CqAq - boAia + Ci = 

Equating expressions in equation H2.12|l at different degrees of Y to zero we 
have algebraic equations for the parameters ao, ai, 60 and Ci in the form 

ao = l (2.13) 



Ai=-^bo (2.14) 
5 

Ao - 4^o' + Co + 8a (2.15) 
25 



a=-bo^ (2.16) 



C.=-^C + ico^ (2.17) 

As a result we obtain the nonlinear second order differential equation in the 
form 



1 18 1 

yzz + ^y' - boy. -Coy-— feo' + ^Co' = (2.18) 

Exact solution of equation H2.18|l is determined by the formula 
y{z) = -^bo^ + Co-^boY-UY^ (2.19) 
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Where Y{z) is a solution of the Riccati equation 



Y. + ^ (2.20) 

which is 



y(z)=±^tanh(^±^ + (^o) (2.21) 

Here (po is arbitrary constant. Actually we find equation H2.18(l with fasci- 
nating history. The matter is this equation can be found from the Korteveg — 
de Vries - Burgers equation which takes the form 



ut + uux + Uxxx = boUxx (2.22) 

Equation 12.22|l is a generalization of the Korteveg - de Vries equation in 
the case of dissipative processes Equation (|2.22|l in contradistinction to 

the Korteveg - de Vries equation is not integrable equation. 

Using the travelling wave H1.2|l we have from (|2.22|l after integration over z 
the following equation 



Ci - Coy + -y" + v.. - h^y. = (2.23) 

We can see that equation H2.23|l is equation H2.18(l as well. 

However exact solutions of equation (|2.23|) were unknown over many years. 
We hope that exact solutions (|2.19() of equation (|2.17() were found first in the 
work jllj . However this exact solution have rediscovered in a number of papers 

EilSSlESI 

Assuming &o = in equation l|2.18|l we have exactly solvable equation so- 
lution. This one can be found via the elliptic function. Nonlinear evolution 
equation (|2.22|l in this case takes the form of the Korteveg - de Vries equation 
and solution (|2.19|l is the soliton. 

One can see we can obtain nonlinear integrable ODEs too using our ap- 
proach. 

We have to note that equation (|2.18|) can be transformed to nonlinear ODE 
that is found from the nonlinear evolution equation 



Ut = Uxx - Cqu - dou^ (2.24) 

This is the Fisher equation or the Kolmogorov - Petrovskii - Piskunov 
equation . Special solutions of this equation were found first in the work jSH] 
and repeated a lot of times later. 
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3 Nonlinear ODEs with exact solutions of the 
first degree singularity 

Taking into account nonlinear ODEs (|2.5I) - H2.8() one can find nonlinear differ- 
ential equations with exact solutions. These exact solutions are expressed via 
the general solution of the Riccati equation 

i; = -y2 + a (3.1) 

First of all let us assume that nonlinear ODEs have special solutions of the 
first degree singularity 

y{z) = Ao + A^Y{z) (3.2) 

where Y{z) is a solution of the equation H3.1|l . 

Without loss of generality we can assume = and Ai ~ 1. 

3.1. Second order ODEs. Let us find the nonlinear second order ODEs 
with exact solutions H3.2|l where Y{z) is a solution of equation H3.1|l . 

Second order ODEs with solutions of the first degree singularity can take 
the form 

Vzz + aoyVz + aij/^ - boy^, + 6iy^ - aCoy + aCi = (3.3) 

Here ao, ai, 5o, 6i, Co and Ci are unknown coefficients of equation H3.3|l . 
Our problem is to find some relations between coefficients of equation l|3.3|) so 
that equation H3.3|l has exact solution (|3.2|l . 

Substitution of H3.2|l and Eqs. H3.1|l and (|2.5|) into equation (|3.3|l lead to 
relation 

(ai - ao + 2)y3 + (6^ + hn)Y^ + a{ao - Cq - 2)Y + aCi - a6o = (3.4) 
From equation (|3.4|l we have 

fli = ao - 2, 
bi = -bo, 

(3.5) 

ao = Co + 2, 
6o = Ci 

As a result we obtain the second order ODEs with exact solutions H3.2|l at 
Ao = and Ai = 1 in the form 

y,, + ((2 + Co)y - Ci) y, + Coy^ - Ciy^ - aCoy + aCi = (3.6) 
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Here a is arbitrary constant. 
At Co = -2 wG get ODE 



Vzz - Ciy, - 2t/3 _ Ciy^ + 2ay + aCi 



(3.7) 



This equation can be found from the modified Korteveg - de vries - Burgers 
equation 



using the traveUing wave H1.2|l . We have found that equation H3.8|l has solution 
(|3.2(l . Exact solutions of equation was found in work |4()j . 

3.2. Third order ODEs. Now let us find the nonlinear third order ordi- 
nary differential equations with exact solutions (|3.2|) . 

Nonlinear third order ordinary differential equations of the polynomial type 
can have exact solution of the first degree singularity if this one takes the form 

aoVzzz + aiyyzz + a2y'^yz + as.yz^ + a.iy'^ + bayzz + hyyz + b2y^+ 



Substitution (|:12() . H3.1(l . H2.5|l and l|2.6(l into equation H3.9|l leads to the 




(3.8) 



(3.9) 



+doyz + diy^ - aC'oy + aCi = 



equation 



(a4 - 02 + as - 6 ao + 2 ai) + (2 feo + &2 - ^i) 



+ 



+ (— + di — do + 2 030? + 2 aia — 8 aaa) Y'^ + 



(3.10) 



+ (2 boa — bia — a Co) Y + a Ci — 2 a^a^ — d^a + a^a^ = 



We have from equation (|3.10(l 



04 = 02 + 6ao — 03 — 2ai 



(3.11) 



Co 



(3.12) 



di = do + a{a2 — 2a^ + 8ao — 2ai) — 



(3.13) 



&i = 26o - Co 







(3.14) 
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a= — (3.15) 

fli — zap 

As a result we obtain the nonlinear third order differential equation in the 
form 



aoVzzz + (aiy + bo) y^z + a^Vz"^ + {a2y^ + 2ybo- yCo + do) yz + 

+ (a2-a3 + 6ao-2ai)/-Co2/^+ (3.16) 

+ {a2a — 2 a^a + do + 8 aoa — 2 aia) y'^ — a Coy + aCi =0 

Solution of equation (|3.1()l) is found by formula l|3.1|) at Aq = and Ai = 1. 
Assuming ai — 0, a2 ~ 0, 03 = 0, 60 = Co/2 in equation (|3.16() we get 



aoyzzz + ^Coyzz + doyz + 6ao2/^ - Coy^ + {8aao + do)y'^ - aCoy + aCi = 

(3.17) 

This equation can be found from the generalized Kuramto - Sivashinsky 
equation [Tn[T7| 



ut + eu^Ux + (iuxx + lUxxx + Suxxxx = (3.18) 

if we look for solution of equation (|3.18|l using the travelling wave H1.2|l . It this 
case we have from equation H3.18(l 



Ci - Coy + + I3yz + lyzz + Syzzz = (3.19) 

One can see that equation (|3.17|) and H3.18|l have exact solutions that are 
expressed via general solution of the Riccati equation 

3.3. Fourth order ODEs. General form of nonlinear ordinary differential 
equation of the fourth order can be written as the following 



aoyzzzz + aiyyzzz + a2y'^yzz + a^yzyzz + a^yyz"^ + a^y^yz + aGy^ + 



+boyzzz + biyyzz + b2y'^yz + bsyz^ + b^y'^ + hoyzz + hiyyz+ (3.20) 



+h2y + doyz + diy ~ Coay + Cia = 

Substituting equation H2.5|l - H2.8|l into equation (|3.2(J() as a result we have 
nonlinear fourth order ordinary differential equation with exact solutions that is 
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determined via general solution of the Riccati equation (|2.8|l . It takes the form 



aoVzzzz + {aiV + bo) Vzzz + {biV + a^Vz + 02^^ + ^o) Vzz + 

+ {a^y + 63) yz'^+ 

+ {do + a^y^ + yCo + ^22/^ - ya^a + 2 yaia - 16 yaoo) yz + 

+ (2 yasa + 2yho) yz + (-a4 + 2 03 + 6 ai - 24 ao - 2 02 + ag) y^+ (3.21) 

+ {-b3-2bi + b2+6bo)y^+ 

+ {—a^a + a^a + 2 — 6 oia — 2 030; + 24 aga + Cq) ?/'^+ 

+ (2 63a — 62a ~ 8 boct + do + 2 bio) y^ — Coa y + C\a — 
Here a is determined by the formula 

ai=0, «2 = |^ (3.22) 
2bo - 03 

From equation H3.22(l one can find some nonlinear ODE which are useful in 
physics. 

Assuming ao = 1, ai = —10, 02 = 0, 03 = —10, 04 = 0, 05 = 6, bo = 
fei = 62 = ^3 = &4 = do = c^i = 5o = 3i = (^2 = in equation (|3.21() we have 
nonlinear ODE in the form 

Vzzzz - lOy^z - lOyy' + Gy^ _ aCoy + aCi = (3.23) 

This is integrable equation and is found from the modified Korteveg - de 
Vries equation of the fifth order 01] 

d 

+ — {uxxxx - lOu'^Uxx - IOmm^ + 6u^) = (3.24) 

Using the travelling wave p.2|l from H3.24|l we have equation (|3.23|l . 

We can see that equation H3.23|l is relative to the class of nonlinear ODEs 
(j3.2HI . Special solution of equation is determined by the formula H3.2|l at Ao — 
0, Ai = 1, Co = 6a and Ci ^ 0. 

Assuming ao = 1, ai — —5, a2 = 5, a3 = —5, = 0, = 1, bo = bi — 
^2 = ^3 = ^4 = rfo = f^i = c?2 = 5o = 31 = in (|3.21|l we have equation in the 
form lIllESl 

Vzzzz - ^V^Vz + ^UzHzz - ^vvl + - aC'oy + aCi = (3.25) 
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(3.27) 



This is integrable equation too. This one can be found from the Fordy 
Gibbons equation 



d 

ut + ^ {uxxxx - hu^Ux + ^UxUxx - 5uul + u^) = (3.26) 

if we look for solution in the travelling wave H1.2|l . 

We get that equation H3.25|l has special solution H3.2(l at Ao = 0, Ai = 
1, Co = a and Ci = 0. 

We can see that a number of integrable equations are abundant in the class 
of nonlinear ordinary differential equations H3.21|l . 

Assuming ai = 0, 02 = 0, 03 = 0, 04 = 0, as = 0, 63 = 0, 62 =0, 61 = 
0, /io = — % + Saao we have 



aoiJzzzz + boVzzz + [ ^^C'o + 8 aoa ] yzz + doyz+ 



+ (—04 + 2 03 + 6 ai — 24 ao — 2 02 + 05) y'^+ 
+ (-63-261 + 62 + 660)2/^+ 
+ {—a^a + 040 + 2 — 6 aia — 2 a^a + 24 aoa + Co) j/^+ 

+ (2 63a — 62^ — 8 boa + do + 2 bia) — Coa y + Cia = 

This equation can be found from the nonlinear evolution equation of fifth 
order that takes the form 

ut + {au + + 71/ + Su"^) Ux + doUxx + f 8aao - j Uxxx+ 

(3.28) 

~\'hoUxxxx + ^0^2:a;a;a;a; — 

Some variants of this equation are used at description of nonlinear waves 

miisiiisiiiT]. 

4 Nonlinear ODEs with exact solutions of the 
second degree singularity 

Let us find nonlinear ordinary differential equations with exact solutions 

y{z) ^ Ao + A^Y + A2Y^ (4.1) 
where Y{z) is a solution of the Riccati equation as well. 
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Second order class with solutions (|4.1|l wc found before in Introduction. Now 
let us find the third order class of nonlinear ODEs with solutions (|4.1I) . 

4-1. Third order ODEs. Nonfinear third order ODEs with solution (|4.1() 
can have the form 



aoUzzz + aiyy^ + feoy^z + hy + diy - Coy + Ci (4.2) 

Substitution (|4.1(l into equation H4.2|l lead to the algebraic equations for the 
coeflticients of equation (|4.2|) . As a result of calculations at A2 = —12, Ai = 
12/3, Ao = we have 



ai = ao (4.3) 
bo = 2bi + h(3ao (4.4) 
do = {0^ + 8a)ao + 10/36i (4.5) 



^ 2(/?^ + 8a) 

Using values for a from (|4.7f) we get 
&1^^=-|^, b'^'-lPao (4.8) 

4^' = 3/3^ao - ^, 4'^=6;3^ao-^ (4.9) 



6W=5/3ao-^, 5^2)=6/3ao (4.10) 



C« = ^/32Co, Cf)= 6/3^^0 + ^ (4.11) 
We have two nonlinear ODEs with special solutions 



aoyzzz + aoyyz + ( 5(3ao - ^ ] yzz 



Coy' , ^„.2 5Co\ r< n 
3/3 ao - — - - Coy + -rr^o = 



6/32 V 3/3 7"-' 2 



(4.12) 



12 



and 



1 2 



(4.13) 



Solutions of equations H4.12|l and (|4.13|l are found by the formula 

y{z) ^ 12(3Y - 12Y^ (4.14) 
where Y{z) is a solution of the Riccati equations 



= -Y^ + (4.15) 
and 

Equations l|4.12(l and (|4.13(l are obtained from nonlinear evolution equation 



Ut + XiUUx + X2UXXX + ^SU-xxxx + Xi{uUx)x + ^bUxx = (4.17) 

if we look for solution of (|4.17() in the form of travelling wave (|1.2ll . 

Equation (|4.17(l was used in of works |48ll49llHn] for description of nonlinear 
waves. Some exact solutions were obtained in papers [JHHSOI- 

4-2. Fourth order ODEs. Let us write the general form of the nonlinear 
four order ODEs with exact solution (|4.1|) . It takes the form 



aoUzzzz + aiyy^z + 02?/^ + a-sy' + 602/222 + hyyz+ 

(4.18) 

+^0^22 + diy"^ + hoyz - aCo + Ci = 

Substitution (|4.1() into equation (|4.12f) at A2 = 1, Ai = 0, = gives 
relations for parameters 

ag = -120ao - laa - 6ai (4.19) 



61 = -1260 (4.20) 



di = 240aao + Saoa - 6do + 8aai (4.21) 
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ho — 8abo 



(4.22) 



do = ITaao + \aa2 - \Co + \oiai (4.23) 
z o 4 



Ci = -18a^ao - a^as + ^a^Co - ^"^ai (4-24) 

We have nonlinear fourth order ODE with exact solution (|4.1I) at Ao = 
0, Ai =Q and A2 = 1 in the form 



aoVzzzz + boVzzz + ( aiy + 17aoQ; + ^020 - ^Co + 1/4 aia ) 2/22 + 



+02^^^ + {biy + 8 5oa) j/, + (-120 ao - 4 03 - 6 ai) (4.25) 
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-Co + —aia + 5 020; + 138 aoa j ~ ^oct y + Ci — 

Values a in equation (|2.8(l is found from equation H4.24|l at given parameters 
Ci, oo, 02 and Co- 

Assuming ao = 1, ai — —20, a2 = —10, 03 = 40, bo — hi — do — di — 
/lo = in equation H4.25|l we have 



Vzzzz - 20yy,, - lOyl + 40y3 - aCoy + Ci = (4.26) 

This is integrable equation. This one is found from the Korteveg - de Vries 
equation of the fifth order j41| . 

d 

wt + ^ {ua^xxx - 20uu^.^ - lOul + AOu^) = (4.27) 

Equation (|4.2(j|) has special solution (|4.1|) at Ao = Ai = 0, A2 = 1 if we take 
Co = 56a and Ci = IGa^. 

Assuming ao = 1, ao = —18, a2 — —9, as — 24, bo = bi = do — di — 
/iQ = in equation H4.25|l we get 



zzzz - l^yyzz + 24y3 _ - aCoy + Ci = (4.28) 

This is integrable equation too. This one is found from the Schwarz - Kaup 
Kupershmidt equation |51II52| 



d 

Ut + ^ {u^xxx - 18MMa:x - + 24u^) = (4.29) 
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Equation (|4.28|) has special solution (|4.1|) at Aq = — 2a/3, Ai = and 
^2 1 if we take Cq = and Ci = -8a^/9. 

We can see again that our class of nonlinear ODEs contains a number of 
integrable equations. 

Assuming ai — 0, a2 ~ and 6i = we have from equation (|4.25(l 



aoVzzzz + boHzzz + [ ^^C'o + ITfloa ) yzz +8boay^ - 120y^ao+ 



(4.30) 

^Co + 138 aoa^ y^ - Coay + Ci = 
This equation can be obtained from nonlinear evolution equation in the form 



ut — 36Qaou^Ua: + ( -Co + 276aao ] uux+ 

(4.31) 

+8aaoUxx + ( ITaag - -Co ) Uxxx + bou^xxx + aoUxxxxx = 



if we look for solution using the travelling wave (|1.2I) . Equation (|4.31|l were 
considered before. Exact solutions of this equation were found in works I11II14I 

E31E1|. 

5 Nonlinear ODEs with solution of the third de- 
gree singularity 

Let us find nonlinear ODEs which have exact solution in the form 

y{z) ^ Ao + AiY + A2Y^ + AsY"" (5.1) 

where Y{z) is a solution of the Riccati equation too. One can see that y{z) 
has the third degree singularity. We have not got any nonlinear second order 
ODE in the polynomial form with exact solutions H5.1|l . However we can find 
the nonlinear third order. 

5.1. Third order ODEs. General form of the nonlinear differential equa- 
tion with exact solution l|5.1|l is enough simple. It takes the form 

aoyzzz + fliy^ + bayzz + doyz - 2Co2/ + Ci = (5.2) 
Using new variables |20j one can write equation 15.2|l in the form 

yzzz + aiy"^ + ayzz +yz'~ 2Coy + Ci = (5.3) 
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Substituting solution (|5.1|l into equation H5.3f) and taking into account equa- 
tions 1)2.81) . H2.5|l and H2.6|l we have 



A3 = 60, ai = 1, (5.4) 



A2 = -ya, (5.5) 



30 15a^ ,^ „x 

Ai = -60a+ (5.6) 

19 152 



We also have two additional equations 



75 n 1200 1965 , 330 1305 . , „ 

-19"^ +^"+ 11552 ^ + 361 - 1444^ 

and 



(5.8) 



cr (o- - 4) (cr + 4) (13 cr^ - 56 + 3040 a) = (5.9) 
From equation (|5.9|l we get 

CTi = 0, (72,3 = ±4, (5.10) 

a=- —a^ (5.11) 

380 3040 ^ ' 

Substituting 1)5.10)1 and 1)5.11)1 into equation ()5.8)) we obtain 

1 11 1 

"1 = "2 = y^, 03,4 = ±^ (5.12) 



16 12 , ^ 

-4,5 ^±-^, -a,7 = ±-= (5.13) 



Using (74,5 and erg,/ we have from equation 1)5.11)) 
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As this take place constant Ci takes 



^(1) _ , 225 (2) _ 2 2475 (3) _ 2 q 
^1 "^° + 6859' ^1 "^°" 6859' "^0-9, 



/-•i.**^ — A n^^'i — 2025 (6) _ „2 



(5.15) 



389017 " ^ 103823 



As result of calculations we find six equations of equation H5.2|l with exact 
solutions. 

Assuming ct = and a = —1/76 we have equation 
225 

y.zz +y^+y.~2yCo + -— + Co^ = (5.16) 

6859 

with exact solution in the form 



45 1 / z 

yiz)-Co + -Y + 60Y^ n^)-±^^tsi±-^^+^o] (5.17) 



At a ~ and a 0:2 = 11/76 we obtain 
2475 

yzzz +y^ + y.-2yCo~-— + Co^ = (5.18) 

6859 

with exact solution 



Assuming a — -i and a = 03 = 1/4 we get 

yzzz + + 4y,, + - 2yCo - 9 + Co' = (5.20) 
with exact solution 



,(z) = ^ + Co-15y-30r' + 60r^ y(z) =±^tanh(±^ + (^o) (5-21) 
At (T = 4 and a = 04 = —1/4 we obtain 



yzzzz + ^yzzz + yz + y^~ 2Coy + - 4 (5.22) 
with exact solution 

y{z) = Co - y + i5y - 3or' + 6or^ y{z) = ±i tg (±1 + <^o) (5.23) 
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Assuming a — 16/a/73 and a — — IG/a/TS we have 



16 2 ^9 2025 

2/-- + "711^- + + 2/' - 2Co2/ + Cl = Q- (5.24) 



(5.25) 



with exact solution 



, , ^ 485V73 1075^^ 120 . 

y{z) = Co H 1 Y —Y^ + 60y-^ 

' " 63948 1168 ^73 



y(z) = tanh {±^j= + ipa 

2^/73 V 2VT3 

At a- = 12/\/47 and a = ae = 12/\/47 we get 



12 900 

+ + + 2/^ - 2C„, + Cl - = (5.26) 

with exact solution 



.(.) = Co + ^ + l^r-^y^ + 6on 

' " 4418 47 47 ' 



Y{z) = tanh {±^-= + .^o 

2\/47 V 2\/47 



(5.27) 



All these equations can be found from the Kuramoto - Sivashinsky equation 
using the travelling wave (fT^ . This equation takes the form |^l56ll57llKHj 



ut + 2aiuux + boUxxx + d^u^x + a^Uxxxx = (5.28) 

Solutions of this equation were first found in work |11| but they were re- 
discovered a few times later. At cr = 4 there is periodic solution of equation 
(lOl [TS|. 

We found that there is only one form of the nonlinear third order differential 
equation with exact solution of the third degree singularity. 

5.2. Fourth order ODEs. General form of nonlinear fourth order ODEs 
with solutions 1)5. 1|) can be presented as the following 



aoVzzzz + aiyj/z + ^o^zzz + biV^ + doVzz + hoVz - Coa y + Ci = (5.29) 
For more simple calculations let us first take = and — 120. We have 

ai — ao (5.30) 
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6i = 



(5.31) 



19 

do = —aoAi + 38aao (5.32) 
60 



19 

ho = ^boAi + 38abo ~ ao^o (5-33) 
60 



Ao^l'-^ + 108^ + ^ + 11^, 6o ^ (5.34) 
5 Oq Oq Oo 3600 Oo 

We also have two values for a 



(1) Ai (1) llAi 

a] = , a, — 

^ 360' 2 1080 

As this takes place we obtain two values of the constant Ci 



360 5o 



(5.35) 



(1) ^ A^ , 

^ 2332800 bo 43740 



= + ^ Ai^feo (5.37) 

1 259200 6o 1620 

Equation with exact solution 

^(^) = -3^^ + ^^^ + ^^«^' (^-^^^ 
takes the form 



(360 aoybo + 76 Ai&o^ + ^laoCo) yz 
aoVzzzz + ooj/zzz + doyzz 



(5.39) 



Where Ci is determined by formulas (|5.36|l and l|5.37|l . 
Assuming 6o = we have Ao ^ Q but we find Co in the form 



_ ao (5040 v4ia + 388800^2 + 11 Ai^) 
° ^ 3600 a 
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In this case we have 



(2) 11^ (2) ^1 (2) ^1 

a] = , Q!n — , a, — 

^ 1080 ^ 360 ^ 12C 

and three values of constant Ci in the form 



(5.41) 



C«=0, cf^^O, C['^ = ^^^ (5.42) 

900 

Equations with exact solution 



y{z)^ A0+A1Y+ 120 (5.43) 
can be presented in the form 

aoVzzzz + aoyyz 2J0~'^^^ ~ ooAoj/^ = (5.44) 

aoVzzzz + aoyyz H ^ — yzz - 00^02/2 = (5.45j 

aoyzzzz + aoyyz - aoAoy^ - -^^y + Ci = (5.46) 

We have not got any information about physical application of equations 
(EH, J05ll and (IQIHl . 

6 Nonlinear ODEs with exact solutions of the 
fourth order singularity 

Let us find nonlinear ordinary differential equations which have exact solutions 
of the fourth order singularity. These solutions can be presented by the formula 

y {z) = Ao + AiY + A2Y^ + AsY'^ + AiV^ (6.1) 

where Y{z) satisfies the Riccati equation again. 

We can not suggest nonlinear ODEs of the second and third order of the 
polynomial form with solution H6.1() . In this case we can take the fourth order 
ODE in the form 



yzzzz + aiy"^ + boy^^^_ + doyzz + eaVz ~ C^y + Ci = (6.2) 

For calculations it is convenient to use first A^ = 0, A^^ — 840/3/11 (/3 is 
new parameter), A/^ — — 840n. Substituting H6.1|l into equation H6.2|l we have 

ai = l 6o = /3, do = Y^/3' + 104a (6.3) 
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en = /3 H /3 a H 

" 1331 11 140 



..oro^ 32760 /J^a 

Ai = 443520 — -;r + 



-30008q! + 65/?2 11 -30008 a + 65/32 

11340 /j^ 
~ 1331 -30008 a + 65/32 

We also have equation for a. in the form 

(-12100 a + 27/3^) (-484a + /32) x 

X (81 /3^ - 696960 /3''a + 289774672 /32a2 + 7029554048 oF) 
Assuming 



we also have the constant C\ in the form 



2 



900 , Co 
^ 214358881 4 

In this case we get equation 



Vzzzz + /3 fez + ^ /3^yzz + i^^Vz + y^- 



900 o Co^ „ 
-^°^- 214358881^ +^=« 

with exact solution 



,(.) = Co + ^/3^-^/33r + ^^y3_84oy4 

''^ " 29282 '^ 133r 11 
where 



Y{z) = ±|rtanh ( ±^ 



22 V 22 ' 
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Assuming Ai — 0, A2 — lAOn, — and A4 — —840 we have after 
calculations 



ai = 1, 5o = 0, do = — 39n + 104 a, eo = 0, 



C Q3 

= y + 816 + y - 528 na 



(6.13) 



We get three values for a 



W31 zV3T\ 

a, = -, a,,3 = -(^-±^J (6.14) 
and the constant Ci in the form 



Ci = -324n4 + ^^ (6.15) 
In this case we obtain equation in the form 



Vzzzz + %/ - 13 ny„, - Coy - 324 + ^ = (6-16) 
with exact solution 



where 



Y{z) = ± ^tanh (^±^ + ^0 j (6.18) 

Exact solutions of equation Hfi.l6|l were found in works |13lll4j and rediscov- 
ered in a number of papers later. 



7 Nonlinear ODEs with exact solution of the 
fifth degree singularity 

Assume that nonlinear ODEs have exact solution of the fifth degree singularity. 
Simplest case of this solution takes the form 

y [z] ^Ao+ AiY + A2Y^ + A^Y^ + A^Y" + A^Y'^ (7.1) 
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One can see that nonlinearity for equation with solution (|7.1|) is . We have 
term of the tenth degree singularity of this equation. If we want to have the 
polynomial form of this equation we have to take 



Uzzzzz+aiy =0 (7.2) 

However we can add to equation (|7.2(l other terms with lesser singularity. 
As a result we have nonlinear ODE in the form 



+ ^o2/zzzz - rfoyzzz + c^y zz + h^y z + aiy"^ ~ Coy + Ci = (7.3) 

This discussion were given by N.A. Kudryashov in the work 26_ to look for 
the truncated expansions for nonlinear ODE with the simplest nonlinearity. 

Equation H7.3|) is found from nonlinear evolution equation which takes the 
form 



'2ciiUUx ~t~ hoUxx ~^ ^O^xxx ^qUxxxx ~t~ ^O^xxxxx ~t~ Uxxxxxx — ('^■'^) 

Last years equation (|7.4|l was used for the description of the chaos model 
|59[l6l)[ICTll62| and it is important to have exact solutions of this equation. Let 
us find exact solutions of equation (|7.3|) at eo = and bo = Q. 

Using new variables 



(7.5) 



C'o = Codo C'l = CicIq ^^^B ^ 
We have equation from H7.3|l in the form 



Vzzzzz - yzzz + <yyz + ^y^ - Coy + Ci = (7.6) 

(primes are omitted in equation (|7.t)|) 'l. 

Substituting H7.1|) into equation H7.6|l we obtain 



2520 

A5 = 30240, A4 = 0, A3 = — - 50400a, A2 = 0, 



1260 2 12600 2520 

^^ = ^^ + ™a^- 30371 +^"' ^° = ^" 



(7.7) 



We also have 



_ -92400 a + 10204656 + 213811840 + 2045 
^ ^ 121(9240a- 79) ^ ' ' 
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Taking into account (|7.8|l we get six values for a 



46031 557 

a4 



52800m 52800 52800' 

TO 46031 557 I to 46031 

"^'^ ~ 105600 ~ 105600TO ~ 52800 "ss" V 1200 ~ 1200to 

Where 

TO = (ll3816753 + 1260V8221079733) ' 
Assuming a = ai we have equation 

3259 1 2 ^ 321489 1^ ^ _ 

with exact solution 

189 3780 

where 



2/ (^) = C^o + ^ 5^ - ^ + 30240 



1 ( z 

Yiz) = ± — ^= tanh ± — r= + osn 

2\/TTo V 2\/TTo 

At a = q;2 we obtain 



1095 1 2 ^ 12403125 1^ 2 _ 
y..... - 2/... - + 2?/ " ^02^ ~ 82458112 + 2^° 

with exact solution 



where 



Assuming a = as we have equation 



114 1 2 127008 2 

fc...-y... + ^2/.+ 2y -2 ^"2^ - 20131375 +^^° ^ 
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with exact solution in the form 



y{^)-C,+ '^Y^ + 30240 r^ (7.17) 



Y{z) = tanh (±^1= + (fo) 

2V55 V 2^/55 J 

We hope these exact solutions will be useful at the study of the turbulence 
processes where equation (17.411 is used. 

8 Conclusion 

Let us emphasize in brief the results of this work. We noted that a lot of 
nonlinear differential equations have exact solutions expressed via hyperbolic 
functions. It is well known that the base of the hyperbolic functions is the 
general solution of the Riccati equation. The widely used tanh method applied 
for finding exact solutions of many nonlinear differential equations confirms this 
idea. 

This observation suggested that a lot of special solutions of nonlinear dif- 
ferential equations can be presented in the form of the general solution of the 
Riccati equation. Using this observation we formulated new problem that is to 
find nonlinear ODEs with exact solutions. We have found a number of nonlin- 
ear differential equations of the second, third and fourth order which have exact 
solutions. Exact solutions of these equations have different singularities and are 
expressed via general solutions of the Riccati equation. We also list a number 
of nonlinear ODEs with exact solutions that are found from the widely used 
nonlinear evolution equations. 
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